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Feynman integrals in the physical light-cone gauge are harder to solve than their covariant 
counterparts. The difficulty is associated with the presence of unphysical singularities due to 
the inherent residual gauge freedom in the intermediate boson propagators constrained within 
this gauge choice. In order to circumvent these non-physical singularities, the headlong 
approach has always been to call for mathematical devices — prescriptions — some successful 
ones and others not so much so. A more elegant approach is to consider the propagator from 
its physical point of view, that is, an object obeying basic principles such as causality. Once 
this fact is realized and carefully taken into account, the crutch of prescriptions can be 
avoided altogether. An alternative third approach, which for practical computations could 
dispense with prescriptions as well as prescinding the necessity of careful stepwise watching 
out of causality would be of great advantage. And this third option is realizable within the 
context of negative dimensions, or as it has been coined, negative dimensional integration 
method (NDIM). 

■ §1. Introduction 

Light-cone gauge in quantum field theory is a fascinating subject. It is a "simple" 
gauge choice in the sense that the emerging propagator for the boson field has a 
Oh, deceivingly simple structure, when compared to other non-covariant gauge choices. 

Yet, as it has soon been realized, it hid subtle complications when one actually 
wanted to work with it. 

To make things more concrete, let us analyse it in the framework of vector gauge 
fields, e.g. the pure Yang-Mills fields, where, after taking the limit of vanishing gauge 
C3 ■ parameter, the propagator reads, 



(N 
> 



k 2 + ie 



n u.ky + n u kn 



k-n 



where is the arbritary and constant null four-vector which defines the gauge, 
n-A a (x) = 0; n 2 = 0. This generates Z)-dimensional Feynman integrals of the 
following generic form, 



Lie 



d D ki f(kj-n*, Pr n*) 



A(kj, pi) h(kj-n, p r n) 
where p\ labels all the external momenta, and n* is a null four- vector, dual to n M . 
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We will evaluate, as a pedagogical example, an integral where, 



A(kj, Pi) = {k 2 )- i (k-p) 



-J 



h(kj-n, pi-n) = (k-n) l , f(kj-n*, pi-n*) = (k-n*) T 

(1-3) 



with negative and m positive or zero. 

A conspicuous feature that we need to note first of all, is that the dual vector 
n*, when it appears at all, it does so always and only in the numerators of the 
integrands. And herein comes the first seemingly "misterious" facet of light-cone 
gauge. How come that from a propagator expression like (1-1), which contains no 
n* factors, can arise integrals of the form (1-2), with prominently seen n* factors? 
Again, this is most easily seen in the framework of definite external vectors n and 
n*. An alternative way of writing the generic form of a light-cone integral is 



1 lc 



J 



dPk; 



Pi 



Ml ^ 



(1-4) 



A(kj, pi) h(kj-n, pi-n)' 

where the numerator g(kj J , pf) defines a tensorial structure in the integral. For a 
vector, we have k^ = (k + , k~ , k T ), where k + = X(k° + k D ^ 1 ) and k~ = X(k° — k ^ 1 ) 
with A being a normalization factor. If we choose definite n and n* such that = 
(1, 0, 1), and n* = (1, 0, —1), this allows us to write k + = k-n and 

k~ = k-n*. We have therefore traced back the origin for the numerator factors 
containing n*. By the way, these terms have nothing whatsoever to do with some 
kind of prescription input as it is in (2-5) below. There, the n* factors are added 
into the structure by hand via the ad hoc prescription. 



§2. Prescriptions 

Now, the troublesome factors in the denominator, represented by h(kj-n, p\-n) 
for one-loop, two-point functions are typically products of the form 



(2-1) 



h(kj -n, pi -n) (k ■ n) [k — p) ■ n 



In the standard approach, these are dealt with first by partial fractioning them, 
a trick sometimes called "decomposition formula" (see, for example, ^ ' 3 )) 



(k-n) (p — k)-n 



1 

p-n 



1 

k-n 



+ 



1 



p-n ^ 



(2-2) 



(p — k)-n_ 

It is easy to see that such kind of trick becomes very clumsy to deal with when 
one has higher powers of denominator factors. In fact, the authors of 3 ) agree that 
the application of (2-2) complicates considerably the calculation of such Feynman 
integrals. For example, consider, 



(k-n) 2 (p — k)-n (p-n)" 2 



p-n 1 1 

— 1 1 

(k-n) 2 k-n (p — k)-n 



p-n^0 (2-3) 
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and so on and so forth. Yet, the standard approach based upon prescriptions, cannot 
handle denominator products like these without first decomposing them. As far 
as we know only Leibbrandt and Nyeo said explicitly before, that "decomposition 
formulas" like (2-2) and (2-3) are in fact part of the prescription 3 ) . Be it as it 
is, after the decomposition is made, then and only then, one seeks to handle the 
isolated denominators calling for prescriptions that would turn them manageable 
mathematically. 

Early efforts in this direction draw heavily from the use of the Cauchy principal 
value (CPV or PV for short), 

PV (—] — ► - lim ( — - + ] —) . (2-4) 

\q-n J 2t^0\q-n + ie q-n — it) 

Unfortunately this prescription, although preserves the general structure of the light- 
cone integral (1-2), and is mathematically consistent, does not provide physically 
acceptable results in the light-cone gauge 4 )' 5 ). 

Later on, independently Mandelstam 6 ) and Leibbrandt 7 ) proposed new pre- 
scriptions to treat the gauge-dependent poles, 

1 1 

q-n q-n + ie sgn(q-n*) 
q-n* 
(q-n) (q-n*) + ie 

where the former is the modified Mandelstam prescription and the latter is the Leib- 
brandt's one. It can be proven that, in fact, they are completely equivalent to each 
other, so that sometimes people refer to this as the Mandelstam-Leibbrandt (ML) 
prescription. A conspicuous feature of these prescriptions is that they introduce the 
factor (q-n*) in the denominator of integrals, thus modifying the original structure. 
This, however, does not seem to affect the end result, as far as physically meaningful 
results are concerned. 

Applying the ML prescription for example, in the one-loop two-point function 
integrals, double poles still appear but they cancel against each other and one is 
left with physical poles only. In other words, those singularities that are merely 
gauge artifacts cancel out when one uses this approach. Experts 5 ) argue, therefore, 
that adopting the ML prescription makes the light-cone gauge acceptable, at least 
perturbatively. Along this line, some two-loop calculations have been performed 
explicitly and can be found in the pertinent literature 8 \ 

Pimentel et al. realized that when the physical principle of causality is correctly 
taken into account — either by careful, stepwise watching out for it to be preserved 
along the whole computation, or by implementing it directly onto the propagator, via 
causal considerations the same picture arises, namely, unphysical, gauge-dependent 
poles cancel out leaving solely physical ones 9 - 1 ' 10 ). Actually, the procedure corre- 
sponds exactly to taking out the zero-frequency mode from the Fourier series expan- 
sion for the field operators, ensuring that only positive energy quanta propagate into 
the future and vice versa. Mathematically this can be acomplished via 



(2-5) 
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<^V= PV (^W) - »^ U - 1) (^>^(A (2-6) 

where PV stands for Cauchy principal value. The second term of this expression is 
crucial. It is the term that guarantees that the troublesome zero-frequency mode of 
the field quanta is subtracted out from the energy spectrum. It is exactly tantamount 
to avoiding the pathological mixing of quanta of opposite energy signs propagating 
into the future, thus violating causality. 

By the very nature in which these prescriptions are introduced, they are com- 
pletely powerless to deal with products of light-cone poles without ere using the 
partial fractioning trick. Moreover, in addition to the need to use the "decomposi- 
tion formula", the necessity to resort to prescriptions or even clever maneuvers in 
order to extract out the culprit factor from the original integral is an awkward trail 
to walk along to reach our objective. Because a prescription is solely a prescription, 
a mere device to by-pass a problem, prescriptions are usually inconvenient and unde- 
sirable. Yet, as awkward and burdensome as it might be, there was no other way we 
could do things properly in dealing with perturbative light-cone gauge computations. 

§3. Negative dimensional approach 

Not so now, with the advent of NDIM 11 ). NDIM is a technique wherein the 
principle of analytic continuation plays a key role. With it we solve a "Feynman-like" 
polynomial fermionic integral, i.e., a loop integral in negative D-dimensional space 
with propagators raised to positive powers in the numerator. Solutions arise as linear 
combinations of solutions for simple systems of linear equations and these are then 
analytically continued to allow for negative values of exponents (i.e. propagators now 
become raised to negative powers, becoming denominator terms in the integrands) 
and positive dimension 12 ) . This procedure, of course, when applied to the light-cone 
case has to consider the general structure (1-2) above which is characteristic of this 
gauge. That is, terms like f(kj-n*, pi-n*) ~ (q-n*) a [(q — p)-n*] b which can appear in 
the numerator, evidently must remain there; therefore, no exponents of such terms 
are to be analytically continued to allow for negative values 13 ). Schematically (see 
for instance second reference in 12 ) ) , 

/ d D h A(kj, n) f(kjn*, p^Mkjn, ptn) ^[ f /j^'"*' ^ (3-1) 
J J A(kj, pi) h(kj-n, pi-n) 

where the left-hand side shows the negative dimensional integral to be performed 
and the right-hand side displays the generic form (1-2) for the light-cone integrals. 
Note that the factor f(kj-n*, pi-n*) remains in the numerator of the integrands in 
both sides. The left-hand side of (3-1), which is evaluated via NDIM methodology, 
is defined from projecting out powers of Gaussian type D-dimensional momentum 
integrals of propagators n ) . It is worth mentioning here that the usual Schwinger 
exponentiation of propagators for the factor [h(kj-n, pi-n)] -1 in positive dimensional 
calculation does not work for light-cone integrals 10 ) . 
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Now, one could rightfully ask: how can one get the standard ML results for one 
and two loop ML dimensionally regulated Feynman-integrals with the procedure of 
negative dimensions? Substitute (1-3) in (1-2), that is, consider the integral, 



B(i,j, l,m)= J d D q (q 2 y [(q - p) 2 \ (q ■ n)\q 



n 



(3-2) 



Let us evaluate it in great detail following the steps thoroughly described in our 
previous papers 12 ). Let our starting point be, 



G 



q exp 



-aq 2 — P(q — p) 2 — jq ■ n — Oq ■ n" 



D/2 



exp 



1 



1 



— ( a/3p + j3^jp ■ n + fiQp ■ n* — - jOn ■ n 



, (3-3) 



where A = a + (3. Taylor expanding the exponentials above and solving for the 
integral B(i,j,l,m) we obtain, 



B(i,j,l,m) = {-TT) D / 2 i\jW.m\(-a - D/2)l ]T 



{X,Y}=0 



{p 2 ) X ^{p + ) X *{p-) X * 

Xi!X2!X 3 !X 4 !li!Y 2 ! 



-nn 



x 4 



xS a i5b,j5 c iSd,m5 e „, 



(3-4) 



where a = i + j + l + m + D/2, a = + Y u b = X± + X 2 + X 3 + Y 2 , c = X 2 + X 4 , 
d = X% + X4 , e = X\ + X 2 + X3 + X4 . This system of linear algebraic equations has 
six possible solutions - since the number of unknowns is bigger than the number of 
equations. In all six cases eliminating the deltas leaves us with one remaining sum, 
a hypergeometric series 3F2. Three out of these six series have as its variable, 



z = 



p 2 nn* 
2p + p~ 



while the other three have as variable 



w = z 



2p + p 
p 2 nn* 



In our previous work 13 ) , we have considered the simplest one-loop integral in the 
light-cone gauge, evaluating it via NDIM method, with one- and two-degree violation 
of Lorentz covariance (more specifically, with n M only and with both and its dual 
n*). For the first case we obtained the usual PV prescription result, whereas in the 
second case we obtained the Mandelstam-Leibbrandt prescription result. However, 
that toy integral is easily evaluated in the NDIM scheme, since there are no left 
over summation indices, that is, the number of linear equations in the system equals 
the number of unknowns and the result can be expressed just in terms of product 
of gamma functions. The integrals we consider here are more complicated; they 
generate several linearly independent and dependent solutions which need to be 
sorted out carefully. That is, the space of solutions spanned by base functions are 
not all linearly independent ones. This is characteristic of these particular cases of 
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integrals we are considering here, and this is the reason we treat them at lenght and 
in depth. 

From our previous works 12 ) we know that the Feynman integral will be rep- 
resented by a linear combination of linearly independent series. Moreover, these 
representations are related through analytic continuation. Let us consider the rep- 
resentation of B(i,j,l,m) in terms of z^ 1 , 

B(i,j,l,m) = B 1 + B 2 , (3-5) 

where 

D/2 r(i + j)r(i + i)r(i + m )(p>y(p+y+i+W(p-y+™+w 

1 1 7T) r(l+j + l + D/2)r(l+j + m + D/2)r(l-j-D/2) 

/_ nn *\-3-D/2 

x ( 2 j 3-^2(01, 6i,ci;ei,/i|ii7), (3-6) 

and 

D/2 r(i + i)r(i + j)r(i + m )r(i -a- p/2)( P 2 r- m ( P -)- l + m 

2 [ n) r(l-j-l-D/2)r(l-i-m-D/2)r(l-l + m)r(l + a-m) 
— nn* \ ' 

— 2~ j 3^2(02, b 2 , c 2 ;e 2 ,f 2 \w), (3-7) 

where ai = -i, 61 = D/2+j, a = a + D/2, e x = l+j+l+D/2, fi = l+j + m+D/2, 
a 2 = -I, b 2 = -i-j-l-D/2, c 2 = i+m+D/2, e 2 = 1-l+m, f 2 = l-j-l-D/2 and 
3F2 is the generalized hypergeometric function 14 ) . This is the result in the negative 
dimension region and i,j,l,m positive. Now, in order to be physically meaningful 
it must be analytically continued to positive dimension and I negative, while 
the exponent m must be kept untouched. To carry out this analytic continuation 
one must rewrite the gamma factors as Pochhammer symbols and use one of its 
properties 

Wto'-^' <■!-'■>= (T^BT (3 ' 8) 

Performing these simple operations one gets, 

B AC (i j I m) = vr^ 2 til ~ 1 ~ D l 2 ^ ± 1 ± D M tfyb+y+i+Dft 
[ ,J ' ' > (l + m\j + D/2) 

x{p - )j+ m + D/2 (^*_y J - D/2 3F2 ( ai)6l)Ci;ei)/lH 

T D/2 {-3\i + 3 + m + D/2)(-i\i + 3 + I + D/2) 
(1 + m\ - I) 

x{a + D/2\ - 2a - D/2 + m)(p 2 ) a - m (p-)- l+m 

x 3 F 2 (a 2 ,b 2 ,c 2 ;e 2 ,f 2 \w). (3-9) 

In fact, from the theory of hypergeometric functions 14 -*, we know that the hy- 
pergeometric differential equation for p F q has up to p linearly independent solutions, 
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so, in writing down B(i,j,l,m) we would expect it to contain up to three terms. 
Indeed, among the three series of funtions 3F 2 (- • • \w) there is a third term, namely, 

D/2 r(i + grp. + j)r(i + pr(i -*- D/2 )^y+^+^ 

3 1 n) r{l + i + j + m + D/2)r{l + l-m)r(l-j-m-D/2) 

f +\l-m / —nn*\ m 

x ni-i-i-D/2) (— ) (340) 

where 03 = —i — j — m — D/2, 63 = —m, C3 = i + I + D/2, es = 1 + I — m, 
fa = 1 — j — m — D/2. However, this solution is not a linearly independent one. 
Indeed, from the fact that m is always positive or zero we can rewrite B3 using 
eq.(6), section 3.2 of ref. 1A > — if we choose a.\ = C3, Q2 = 03, n = — e%, p\ = fa - 
as a linear combination of B\ and £?2- 

The integral we are considering here was studied by Lee and Milgram 15 ). Our 
result can also be written in terms of the so-called G— function or Meijer's function 
and agrees with their calculation. Taking the special case where i = j = I = — 1 and 
m = one gets, 

n( 1 1 1 m ( -^d/2 ^(2 ~ D/2)r(D/2 - 1) ( p +)g/2-2 (p - ) g/2-i 
-1,-1,0) = -(-*) - 2 

9 \ 0/2-1 

- J 2^1(1, D - 3; D/2\w) - n D / 2 r(D/2 - 1) 



r(D/2 - 2)r(3 - D/2) ^ D/2 _ 3 ^_ ( -2 



r(D - 3) yiJ ' P \nn* 

x 2 F 1 (l,D/2-l;2\w). (3-11) 

Observe that when D = 4 — 2e there are simple poles in these two terms but they 
cancel since, 

hmr(e) + r(- e ) = 0(l). 

To extract the finite part one can proceed in the same way as we did in the second 
paper of ref. 12 \ The other representation for the integral B(i, j,l,m) also coincide 
with that presented in 15 ). In fact, they are related through analytic continuation (see 
ref. 14 ) for such formulas) - as all representations provided by NDIM for Feynman 
integrals in covariant 12 ) and non-covariant gauges 13 ) . 

Consider now an integral containing two factors of the form (k ■ n), 



T 2 (i,j,l,m) = J d D q(q-p) 2i (q ■ n)' [(q - p) ■ nf (q ■ n*) m , (3-12) 

where the exponent m > 0. The standard approach makes use of decomposition 
formulas (2-2) and (2-3). On the other hand, NDIM can solve all integrals of such 
kind at the same time, giving in positive dimension 17 \ 



D/2+i 

T 2 (i,j,l,m) = 7r 



D/2 1 2p+p- \ r{i + l + D/2)r(l + m) 



n-n* I (p + )- j - l (p~)- m r(-i)r(-j) 
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X jr( l + i + m + jD/ 2)' V L6 > 

observe that the Pochhammer symbol which contain (1 + m) was not analytically 
continued since m must be either positive or zero. 



§4. Conclusion 



Finally, in concluding this work we would like to emphasize some important 
points concerning the application of NDIM technique to light-cone integrals. First 
of all, why one needs to define the original Gaussian integral (3-3) with two-degree 
violation of Lorentz covariance in order to get "causality" preserving results? The 
answer is related to the very definition of light-cone gauge condition where one 
has the external four- vector n^ such that it is chosen to be light-like, i.e., n 2 = 0. 
However, the "light-likeness" condition n 2 = does not uniquely define the needed 
external vector n^ to implement the gauge condition. The reason for this is more 
easily seen considering a particular case where n M = (no, 0, 0), in which case the 
condition n 2 = gives as solutions either no = +n3, or no = — with no > 0. 
Therefore, the components of the "light-like" vector are not linearly independent; 
hence the two possibilities: either n M = (no, +n^, 0, 0) or n* = (no, — n%, 0, 0). 
These peculiar light-cone properties have been shown by G. Leibbrandt 16 ) to be 
connected to the Newman-Penrose tetrad formalism in the context of gravitation 
and cosmology, where a four-dimensional basis is spanned entirely by null vectors. 
In his work, Leibbrandt demonstrated that the two-dimensional vector sub-space 
(no, nz) cannot be spanned solely by the vector n M ; it is not sufficient for this purpose 
because this vector possesses linearly dependent components. A thorough discussion, 
including not only the sub-space (no, n^), but the entire four-dimensional space, is 
found in the reference above given. 

This is the reason why our NDIM calculation done in 13 ) with only one degree 
violation of Lorentz covariance failed to be "causal", reproducing the well-known 
"causality" breaking PV result. The cure for this pathological result could only be 
achieved by introducing the needed additional n* vector, so that the basis for span- 
ning the entire four-dimensional vector space be unique and well-defined. Therefore 
our calculations here always included this needed dual vector n* , so that our results 
are well-defined and unambiguous. 

This NDIM technique is therefore the most powerful and beautiful machinery 
ever to come to front until now to handle the light-cone Feynman integrals. All 
the desirable features are embedded in it and we can list them as follows: i) No 
"decomposition formula" like (2-2) is needed to separate gauge dependent poles 
before the integral is evaluated; ii) It preserves the general structure of the light-cone 
integral, namely, no factors containing n* are introduced in the denominators; Hi) No 
prescription is needed to handle {q ■ n)~ a = type singularities to solve the integral; 
iv) There are no parametric integrals to perform; v) There is no need to split the 
dimensionality of space-time to work out integrations over component sub-spaces like 
in 2 ); vi) Results are obtained for arbitrary negative exponents of propagators and 
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gauge-dependent poles, so that special cases - which agree 1 )' 4 )' 13 - )> with the ones 
calculated using ML prescription - are contained in them; vii) Results are always 
within the context of dimensional regularization, i.e., preserving gauge symmetry. 

The power of NDIM is therefore readily apparent. If one remembers how inte- 
grals over Grassmann variables are introduced, he/she can recall that the underlying 
property employed to define them is translation invariance. It is an outstanding thing 
that this sole property in the fermionic integration is able to guarantee a prescrip- 
tionless light-cone and more important than that, a causality preserving method of 
direct computation of Feynman integrals. 

Of course, like in any gauge, two-loop integrals are more demanding when com- 
pared to the one-loop case, and specially in the light-cone gauge. Preliminary results 
we have show that NDIM can also handle them with more easiness than in the usual 
standard approaches. 
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